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Abstract 

A model is presented which demonstrates that the at- 
tosecond pulse structure of a High Harmonic Generation 
(HHG) seed may be retained through to saturation in an 
FEL amplifier. At wavelengths of ~ 12 nm, a train of at- 
tosecond pulses of widths ~ 300 attoseconds with peak 
powers in excess of 1 GW are predicted from full 3D sim- 
ulation. Methods for improving these results are discussed. 

INTRODUCTION 

High-gain free-electron lasers (FELs) operating at short 
wavelengths have the potential to gather incredibly detailed 
information on how matter interacts and arranges itself at 
the atomic and molecular scales Q. If high brightness 
pulses of sufficiently short duration (10~^^ — 10~^''s) and 
power can be generated, it may also be possible to capture 
ephemeral events, such as molecular bond formation, with- 
out temporal blurring. While HHG delivers trains of at- 
tosecond pulses, they currently lack the peak powers avail- 
able from FELs. 

Several techniques to achieve such pulse durations have 
been proposed (see e.g. lEl [3l IH and references therein) 
where the widths of the pulses generated are of the order 
of the coherence length Ic = A^/47rp, where is the res- 
onant wavelength and p is the fundamental FEL parame- 
ter 0. 

A different technique, that applies the concepts from 
mode-locked cavity lasers, has also been proposed |6| that 
may generate a train of pulses with lengths less than Ic. In 
this technique, a series of spatiotemporal shifts are intro- 
duced between the radiation and the co-propagating elec- 
tron bunch that define a set of axial radiation modes which 
couple via the FEL interaction. The spatiotemporal shifts 
are achieved by delaying the electron bunch using mag- 
netic chicanes inserted between undulator modules. By in- 
troducing a modulation at the mode spacing, each mode 
develops sidebands that can overlap neighbouring modes, 
allowing mode-locking to occur. This is analogous with 
mode-locking in conventional cavity lasers. 

The spectrum generated by the XUV mode-coupled 
SASE system of 1 6 1 has a modal structure with modespac- 
ing Aij similar to the spacing of harmonics in a HHG 
source. Equivalently, the temporal structure of the two sys- 
tems are a comb of attosecond pulses. By matching the 
spectral/temporal structures of a mode-coupled FEL am- 
plifier to an HHG seed, it may therefore be possible to am- 
plify the HHG seed to saturation while retaining its attosec- 
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Figure 1: Schematic of the mode-coupled amplifier FEL 
interaction operating with an HHG seed. The inset shows 
detail of the electron delay. 



ond train structure. Furthermore, it should be possible to 
achieve this without the need for introducing a modulation 
at the mode spacing as the HHG attosecond seed pulses al- 
ready have a well-defined phase relationship. In this paper 
it is demonstrated via ID and 3D numerical simulations 
that this is the case. 

MODE GENERATION IN AN AMPLIFIER 

The spatio-temporal shifts that generate the equally 
spaced modal structure in the amplifier spectrum are 
achieved by periodically delaying the electron bunch us- 
ing magnetic chicanes between undulator modules |6|, as 
shown schematically in Fig. [T] The discrete modal spec- 
trum is similar to the axial modes of a conventional laser. 
These modes are intrinsically coupled via their collective 
interaction with the electron bunch over regions of a co- 
herence length Icoh ~ Ic where Ic = Xr/^7rp is the FEL 
cooperation length. 

The notation used here follows that of ["61 . The total slip- 
page of the radiation with respect to the electron bunch per 
module in units of the scaled electron bunch co-ordinate 
zi |6| is 5 = I -\- S, where I is the slippage occurring 
in the undulator and S is the slippage within the chicane. 
The scaled resonant FEL frequency is O = {co — u;r)/2pu;r 
where ujr is the resonant FEL frequency. The slippage en- 
hancement factor is defined sls Se = s/l. In [6|, a solution 
is obtained for the one dimensional wave equation describ- 



ing the field evolution of a small, constant electron source 
term for a series of N undulator-chicane modules. For 
> 1 the spectrum takes the form of a sinc-function enve- 
lope (the single undulator module spectrum) modulated by 
a frequency comb centered at the scaled resonant frequency 
O, with mode separation AO = 27t/s corresponding to a 
Auj = 27r/Ts where Tg = s/c = slc/cis the time taken 
for radiation to travel the slippage length. The number of 
modes under the central peak of the spectrum is shown to 
be A^o = 2Se-l 

AMPLIFICATION OF AN HHG SEED 

The temporal profile of the HHG seed is a comb of at- 
tosecond pulses separated by half the wavelength of the 
drive laser, A^^. The spectral and temporal structure of 
the radiation generated by the undulator-chicane system is 
matched to that of the HHG seed by setting s = Xd/2, 
where A^^ = Xd/h is the scaled drive laser wavelength. 

For a typical HHG drive laser wavelength of 
Xd = 805 nm (e.g. Ti: Sapphire) and an FEL parame- 
ter of p = 2 x 10~^ similar to that for an FEL operating 
in the XUV, then for undulator modules of 8 periods 
I = 0.201. The value of the scaled slippage generated by 
the chicanes is chosen to be J = 0.616 so that s = Ac^/2, 
matching the undulator-chicane system to the spectral 
and temporal structure HHG seed. This gives a slippage 
enhancement factor of = 4.0625. 

The HHG seed was modelled in a similar way to that de- 
scribed in 171 with the resonant wavelength chosen to be 
the 65th harmonic of the drive laser i.e. A^ = Xd/Q^. The 
Nyquist theorem limits the frequencies that can be mod- 
elled in the averaged code | 7 | to harmonics 33 to 97. The 
temporal and modal structure of the seed is shown in Fig.O 

The system is modelled using a ID code similar to that 
described in (H. The HHG seed is injected into the am- 
plifying undulator-chicane modules with a cold electron 
beam. Note that no energy modulation of the electron beam 
was used. This is different to the mode-locked SASE case 
of f6l which starts up from noise and requires a beam en- 
ergy modulation at the mode spacing Au to couple the 
modes and generate locking between each mode. Here, the 
HHG seed has modes which already have a well-defined 
phase relationship. It will be seen that while this phase re- 
lationship between modes is not significantly altered during 
amplification (the temporal pulse train structure is retained) 
it appears that the lack of modulation at the mode spacing 
Aci; does not allow the modes to 'lock' and the number of 
modes receiving significant amplification is reduced. This 
leads to a broadening of the individual pulses and appears 
in agreement with the scaling observed in 1 6 1 for the pulse 
width Tp (X 1/ ^/Nq where A^o is the number of amplified 
modes. 

After 1 8 undulator-chicane modules the amplified scaled 
radiation power and spectrum is as shown in Fig. [51 It is 
seen that the HHG seed is amplified and the comb of tem- 
poral spikes is retained, but with increased pulse widths as 



discussed. This number of modules corresponds to a FEL 
interaction length of22xr?^4.4. This is just prior to satu- 
ration where scaled powers of | ^ 1 are achieved. How- 
ever, the temporal structure of the pulses begin to break up 
at saturation. (Note that the dispersive strength of the chi- 
canes D = R^q/21c ~ 10 J/6 reduces the gain length as 
discussed in \i6i). 
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Figure 2: Scaled longitudinal intensity profile (top and 
middle) and scaled spectral power distribution for the HHG 
seed (bottom). 



OPTIMISING THE SYSTEM 

In |6| it is demonstrated that the pulse widths in mode- 
locked SASE simulations are proportional to 1 / >/]Vo , this 
is in agreement with an analogous relation from mode- 
locked conventional cavity lasers. 

In order to attain narrower amplified HHG pulses, the 
number of modes under the central peak of the spectrum is 
increased by increasing the value of the slippage enhance- 
ment factor. To increase while continuing to satisfy the 
condition s = Arf/2, the value of Tmust be decreased while 
S is increased to keep s constant. The system is simulated 
for a case with the slippage enhancement Se = 8.125, with 
the number of undulator periods per module equal four. 

The spectral power distributions for the two cases Se = 
4.0625 and Se = 8.125 are shown in Figs.[3]and |4]respec- 
tively. It is seen that the number of modes, A'o, under the 
central peak of the spectrum increases with Se, narrowing 
the pulses. For = 8.125, 34 undulator-chicane modules 
are required to produce the same amplification of the 22 
module = 4.0625 case. 

For larger values of Se, the broadening of the modal 
sinc-function envelope becomes restricted by the Nyquist 
limit. 
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Figure 3: Scaled longitudinal intensity profile (top and 
middle) and scaled spectral power distribution of the am- 
pHfied HHG radiation for = 4.0625, after 22 undulator- 
chicane modules. 

ENERGY SPREAD 

The effects of energy spread on the interaction have been 
investigated using the ID model by introducing a relative 
gaussian energy spread of width cr^/j = p/10 = 0.02% 
for the case Se = 8.125. The results close to saturation are 
plotted in Fig. [51 which is the equivalent of Fig. (H While 
the cold beam case took 34 undulator-chicane modules, the 
case here took 42 modules. Note that the temporal and 
spectral structure of the radiation power remains largely 
unaffected by the energy spread. Larger energy spreads 
have also been simulated with the main effect being a larger 
number of modules to saturation. 

RESULTS OF 3D SIMULATIONS 

The system has also been modelled using the three- 
dimensional code Genesis 1.3 |9|. The main results of the 
ID simulations are reproduced using the 3D code and show 
amplification of an HHG seed pulse retaining its pulse train 
structure. 

Similar parameters as for the XUV mode-coupled FEL 
amplifier of |6| were used (i.e. without electron beam 
modulation) including an energy spread of width a^/j = 
0.01%. The seed power and spectrum show the attosec- 
ond pulse train structure and frequency comb, as seen in 
Fig. [6l The peak power of the seed is 3 MW with FWHM 
duration 33 fs. It was assumed that the wavelength of the 
Ti:Sapphire drive laser was A^^ = 780 nm and the FEL was 
then tuned to be resonant with the 63rd harmonic of the 
drive laser with = 12.4 nm. 
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Figure 4: Scaled longitudinal intensity profile (top and 
middle) and scaled spectral power distribution of the am- 
plified HHG radiation for = 8.125, after 34 undulator- 
chicane modules. 

Each undulator module had 4 periods of 3.1 cm and 
each chicane had a delay of 28 radiation wavelengths cor- 
responding to a slippage enhancement = 8. These pa- 
rameters for the undulator-chicane modules then match the 
mode spacing of the HHG seed. The amplified seed ap- 
proaching saturation, after 48 modules, is plotted in Fig. [7] 
and the pulse temporal and spectral properties agree well 
with the equivalent one-dimensional Se = 8.125 case of 
Fig. [5j The FWHM length of each pulse in the amplified 
train is 300 as. Further 3D simulations (not shown) with 
5'e = 4 agreed well with the Se = 4.0625 case of Fig. [3 

CONCLUSIONS 

It is shown, using ID and 3D simulations, that the at- 
tosecond structure of the HHG seed can be amplified to 
saturation in the mode-coupled amplifier FEL of |6|. It is 
stressed that no pre-conditioning of the electron beam (e.g. 
energy modulation) is required. The only requirement is 
for an undulator-chicane modular structure, the radiation 
spectrum of which matches an HHG seed source. 

The introduction of an energy modulation to the system 
may well improve the amplification of the HHG to generate 
trains of narrower pulses by locking the modal structure. 
This will be investigated in future work. 

Also to be investigated is the behaviour of the system 
when modelled in a non-averaged code, such as ifTOt . Such 
codes have a larger Nyquist frequency and are capable of 
simulating a wider frequency range and so shorter pulse 
lengths. 
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Figure 5: Scaled longitudinal intensity profile (top and 
middle) and scaled spectral power distribution of the am- 
plified HHG radiation for = 8.125 and an initial energy 
spread a^/j = 0.02 %, after 42 undulator-chicane mod- 
ules. 
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Figure 6: Longitudinal intensity profile (top) and spectral 
power distribution (bottom) of the HHG seed in 3D simu- 
lations. 
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Figure 7: Longitudinal intensity profile (top and middle) 
and spectral power distribution (bottom) of the amplified 
HHG radiation in 3D simulations, with = 8. The agree- 
ment with the equivalent ID simulations shown in Fig.[5]is 
very good. 
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